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We study the stability of new neutral and electrically charged four-dimensional black hole solutions 
of Einstein's equations with a positive cosmological constant and conformally coupled scalar field 
The neutral black holes are always unstable. The charged black holes are also shown analytically to 
be unstable for the vast majority of the parameter space of solutions, and we argue using numerical 
techniques that the configurations corresponding to the remainder of the parameter space are also 
unstable. 
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I. INTRODUCTION 

Black hole solutions of the four-dimensional Einstein- 
scalar field system have been extensively studied for over 
thirty years, with a particular focus on proving unique- 
ness ("no-hair") theorems in various models (see, for ex- 
ample, 2] for a review). Within these models, conformal 
coupling of the scalar field is of particular interest, both 
with and without an additional coupling to the Maxwell 
field (see, for example, 3] for a brief review of work on 
this situation). In asymptotically flat space, with no 
scalar self-interaction potential, there is an exact, closed 
form solution (the BBMB solution) 0, H, 01 > which has 
not been without controversy f] because the scalar field 
diverges on the event horizon. Furthermore, it is known 
that this solution is unstable Q. For space-times which 
are asymptotically anti-de Sitter, with zero potential or a 
quadratic scalar field potential, numerical solutions exist 
in four dimensions Q , of which at least some are linearly 
stable under spherically-symmetric perturbations. Inter- 
estingly, the corresponding three-dimensional black hole 
solution is known in closed form 0] , but is unstable ■ 

Asymptotically de Sitter geometries are the focus of 
this paper. For minimally coupled scalar fields, non- 
trivial scalar field hair is possible if the scalar field poten- 
tial is non-convex [ill lia | , although the hair is unstable 
[T^ . For conformally coupled scalar fields, if the scalar 
field potential is zero or quadratic, then there are no non- 
trivial black hole solutions 3J. However, in the presence 
of a quartic self-interaction potential there is an exact, 
closed-form solution for both charged and neutral black 
holes found recently by Martinez, Troncoso and Zanelli 
0, which we shall refer to hereafter as the MTZ solution. 
This solution is the de Sitter analogue of the BBMB so- 
lution, although the scalar field is regular on and outside 
the event horizon. The purpose of this paper is to dis- 
cover whether, hke the BBMB solution, the MTZ black 
hole is unstable. 



The outline of this paper is as follows. In sectionlTTlwe 
briefly review the neutral and electrically charged MTZ 
solutions P], whose stability is then studied using linear 
perturbation theory in sections IIIII and IIVI respectively. 
In the appendix we outline the proof of a result needed 
for the stability analysis of the charged solutions. Finally, 
we present our conclusions in section The metric has 
signature ( — h++) and we use units in which c = SvrG = 
1 throughout. 



II. THE MTZ SOLUTION 

We begin with the action for gravity with a confor- 
mally coupled scalar field with a quartic self-interaction 
potential and an electromagnetic field 



S = 



i? - 2A - g'"'da(l}d^(l) - -R(t)'^ 
6 



(1) 



where a is the coupling constant. For neutral black holes, 
the electromagnetic field is absent. The Einstein and 
scalar field equations derived from this action are 



0; 



(2a) 



□0 - -Rcf) - iact)^ = 0; (2b) 
6 

where □ = g'^^'V^Vj^, and the scalar and electromagnetic 
stress energy tensors are 

+^ [g^^D - + G^„] 0^ - ag^^cj)^; 



T, 



EM 



^ U^'^F^aF.^ - \g^.F^0F'^^ 



For charged black holes, we also have the Maxwell equa- 
tions 
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V^F^'' = 0; 



(3) 
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For both the charged and electrically neutral cases, the 
field equations are solved by the spherically symmetric 
metric 



ds^ = -N{r)dt^ + N{r)-^dr^ + ; 



sin^ 9 dip^ 



where 



(4) 



N{r) 



A 



1 - 



M 



3 \ r 
and the scalar curvature is a constant given by 
i? = 4A. 

The geometry (0J is that of the Reissner-Nordstrom-de 
Sitter black hole with inner, event and cosmological hori- 
zons at values of the radial co-ordinate r given by, respec- 
tively, 



-1 



4M 



1 - 



4Af 



1 - 



4M 
T 



(5) 



where I = y^S/A. From ((SJ, it is clear that the solution 
is defined only for < M < Mmax = 1/4:. 

The form of the scalar field is different for the elec- 
trically neutral and charged models. In the case with 
no electromagnetic field, there is a solution only if the 
coupling constant a is given by 

a = A; 

36 ' 

and then the scalar field takes the form 



(r) 



(6) 



r-M' 

which has a pole at r = M < r+, lying inside the event 
horizon. 

For charged black holes, the only non- vanishing com- 
ponent of the electromagnetic field is 

g. 

2 ' 



Ft, 



-drAt 



where the charge-to-mass ratio is given by 



M / V 36a 



and the scalar field in this case is 



[r = A - 



A M 

6a r — M 



(7) 



(8) 



This latter solution only exists provided a satisfies the 
bound 

36a < -A. 



III. INSTABILITY OF THE MTZ SOLUTION: 
NEUTRAL CASE 

We now analyze the stability of the MTZ black holes 
using linear perturbation theory. The electrically neutral 
and charged solutions need to be considered separately, 
and we begin with the neutral case since this is the sim- 
pler. In this and the following section, we consider only 
spherically symmetric perturbations since these are suf- 
ficient to show instability. Many calculations were per- 
formed using the GRTcnsorll 13] routine linpert 14]. 

The perturbed spherically symmetric metric takes the 
form: 



ds^ - 



^N{t,r)e^^^^*'''Ut^ + N~^{t, r)dr^ 
+r^{d9^ + sm'^ed<fi^); (9) 



where 



N{t,r) = N{r)+eN{t,r); 



(10) 



and e is a small parameter. The scalar field is perturbed 
as 



(j){t,r) = (j){r)+e(j){t,r). 



(11) 



The ansatz (PTT|l can then be substituted into the Ein- 
stein and scalar field equations (|2a|l and (j2b|) . the equa- 
tions then linearized and the metric perturbations N and 
6 eliminated to yield a single perturbation equation for 
the perturbed scalar field cj). While this is possible in 
principle, in practice the algebra becomes somewhat un- 
wieldy, even using a computer algebra package. 

We therefore employ a transformation to a simpler sys- 
tem which makes the calculations more tractable, and 
then substitute back the original perturbations. Under 
the following conformal transformation (see also 
for more details). 



where 



— o2 



= ( 1 - -02 



(12) 



the action (QJ becomes that of a minimally coupled scalar 
field (f> with potential V"(<1>), with 



and 



$ = \/6tanh'^ 



V(^) = 2Asinh2 



V6 
$ 

71 



(13) 



At this stage we should comment that, from Eq. 10, 
there is always a point (r = 2M) between the event and 



cosmological horizons at which 1 - 



0, and therefore 



3 



the conformal transformation H12() breaks down at this 
point. However, since this is only a single point, and not 
an open set, this does not render the theory completely 
ill-defined (see, for example, the discussion in Ref. ^T^). 
In addition, the transformed variables have no physical 
significance; the only purpose of the transformation is to 
simplify the algebra, in particular to derive the master 
perturbation equation, as in the analysis of Ref. 8] . Our 
conclusions will be drawn only once the perturbations of 
the original system have been substituted. 
The new metric H12f) takes the form 

ds^ = -JVe^'^df+JV-^dx^ + x^{d9^ + sm^0dip^); (14) 

where we have transformed the radial coordinate by 

. 1 

x=(l-i^ 

and the metric quantities Af and A are related to the 
original metric variables by Q : 



Af ^ 1- - 



1 ,o 1 



6"^ 6 



'4*'- i'-l"' 



Note that A 7^ in the new metric. The transformed 
variables Af, A and <i> are perturbed as follows: 

Af{t,x) ^ Af{x)+eAf{t,x); 
A{t,x) = A{x) +eA{t,x); 
^{t,x) ^ ^{x) + e^{t,x). 

These expressions are substituted into the new Einstein 
and scalar field equations (a tilde denotes quantities cal- 
culated using the transformed metric (|14|l ): 



^ G^^+A.g,,^- V,,$V^$ + -5^^(V$)2 



= □$ 



dV 



only terms linear in the parameter e are retained, and 
the metric perturbations can then be eliminated. A sin- 
gle perturbation equation is then found for the variable 
(j). For periodic perturbations (f>{t,x) = e^°''^ (f>{x) , this 
equation takes the form: 

-0-Qi^-Q20 + Zi<^ = a2<^; (15) 

where we have introduced the usual "tortoise" co- 
ordinate defined by 



1 



dx Afe^ ' 



(16) 



so that 



dx^ 
dr 



1 



(17) 



and the region between the event and cosmological hori- 
zons r+ < r < r++ becomes —ooKx^, < 00. The 
quantities Qi and Q2 are given by 



Qi 
Q2 



2N 

r 
N 



dN 



NA- 



dV 



dr 2 dr dr 

and the perturbation potential hi takes the form 



U ^ — 



1 - NA-'^B'^ 



dr 



B- 



-Ar' 1 



(18) 



where 



A 
B 
V 



—ref)—: 
6 dr 

1 d^ 
—rip—. 
6 dr 



(19) 



It should be stressed that the equation H15() could equally 
well have been derived directly, without using the confor- 
mal transformation (|12|l . but that our method has simpli- 
fied the algebra considerably. Furthermore, the equation 
(|15|l holds everywhere, although the conformal transfor- 
mation (|12|l is defined only for those values of r for which 
02 < 6. 

The perturbation equation (|15|) has a regular singular 
point at r = ro = 2M, where 0^ = 6. Using the standard 
Frobenius method, near this point cf) behaves like 



<^ - (r - ro) [Ci -I- C2 fog (r - tq)] 



(20) 



for constants Ci and C2 ■ The presence of the regular sin- 
gular point means that it is not possible to define a non- 
singular transformation of l|15|) to standard Schrodinger 
form, as, for example, was possible for three-dimensional 
black holes in Ref. [lOl |. 

In order to cast the perturbation equation (|15|l into 
standard Schrodinger fo rm, we define a new perturbation 
variable ^ as follows [8l[lo|: 



1 



so that the equation (|15|l then takes the form 
52^ 



(21) 



(22) 
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Potential for a black hole configuration inside 9? 




'* Radius 

FIG. 1: Form of the potential U for the neutral MTZ black 
holes. The potential has a double pole at r = ro = 2M. This 
is also the form of the potential for the charged MTZ black 
holes configurations corresponding to points within the region 
K, see section HVI 



with potential U given by H18|l . It is straightforward to 
show that the quantity B H19|l is regular and non-zero 
everywhere between the event and cosmological horizons. 
Therefore, the potential (|18|l is regular for all apart 
from the second term, which has a double pole at r = 
ro = 2M (where 0^ = 6). Near r ^ tq = 2M, the 
potential has the behaviour 

^ = -\ i^iro)? {r - ro)"' + O (r - ro)"' ; 

and — !■ as a;* — > ±oo. The form of the potential is 
sketched in Fig. Like the original equation (|15|l . the 
new perturbation equation H22(l has a regular singular 
point at r = ro, and a Frobenius expansion about this 
point gives the behaviour of ^E* to be 



(r-ro)^ Ci + (72 log (r - ro) 



where Ci and C2 are constants. Using equation (|20|l . it 
can be seen that the transformation H21I) yields precisely 
the correct behaviour for 5*, so that the two perturbation 
equations l(T3|l and are equivalent, in particular they 
have the same number of negative eigenvalues. 

Using a standard result in quantum mechanics ^Sj , po- 
tentials of this form have an infinite number of bound 
states corresponding to < 0, that is, an infinite num- 
ber of unstable modes. Therefore we also have an infinite 
number of unstable modes of the original perturbations. 
This is the same behaviour as found for the BBMB black 
hole l. 



IV. INSTABILITY OF THE MTZ SOLUTION: 
CHARGED CASE 

The perturbation analysis of the electrically charged 
black holes is considerably more complicated, and al- 



though we shall again make use of the conformal transfor- 
mation 1)12(1 to simplify the algebra, the untransformed 
perturbations themselves can be simplified first. 

Making a choice of gauge to set Ar — 0, the gauge 
potential giving rise to the electromagnetic field is per- 
turbed as below: 



At 
Ae 



Q 



+ eAt{t,r)- 



0; 

eA0{t,r); 
eA^{t, r) sa\9] 



so that the second of the Maxwell equations Q is auto- 
matically satisfied. The metric and scalar field are per- 
turbed in the same way as for the neutral black holes 
P jTT|) . The [to), [rO), {tip) and {rip) components of the 
Einstein field equation H2a|) give 



dAeit.r) _ dAe{t,r) _ dA^{t,r) _ dA^{t,r) 



dr 



dt 



dr 



dt 



0; 



and since we can remove an arbitrary constant from the 
gauge potential without changing the resulting electro- 
magnetic field we can, without loss of generality, assume 
that 

Ae = A^ = 0. 

The t and r components of the first Maxwell equation 
© then take the form: 

= yit^r)Q + r^^^Mt^r)+2r^^Mt,r), 



= 



dtdr 



Mt,ry + -^^d{t,r)Q 



These two equations are compatible and can be inte- 
grated to give 



Sit,r) 



Qdr 



where we have absorbed the arbitrary constant of inte- 
gration into the time co-ordinate. With the form of the 
electromagnetic potential now fixed, the {tr) component 
of the Einstein equations can readily be integrated to 
yield 



N = r 



x(l-i, 



1 



(23) 



where T is an arbitrary function of r. Incidentally, equa- 
tion (|23() also holds for the neutral black holes. 

At this point, the perturbation equations are again too 
complicated to readily be simplified and we employ the 
conformal transformation as in the previous section. As 



5 



in section fnH we regard the transformation as a conve- 
nient algebraic tool with no physical significance 0- 

The transformations of the metric and scalar field are 
the same as before p2H13|l , but now the electromagnetic 
field also transforms ilTj] : 



Relations H24I27|) are sufficient to eliminate the met- 
ric and electromagnetic perturbations from the linearized 
scalar field equation. Using the "tortoise" co-ordinate 
(|16II17I) . and following the method of section ITTll we fi- 
nally arrive at the Schrodinger-like equation 



leaving the electromagnetic stress-energy tensor invari- 
ant. The transformed metric variables and scalar field 
are perturbed in the same way as in the previous sec- 
tion. The only non-zero component of the transformed 
electromagnetic field is Ftx, which we write as 



Ftx{t,x)^F{x)+eF{t,x); 



where 



F{x) 



Q 



The electromagnetic and Einstein equations can be 
used to find the following simple expressions for the per- 
turbations of the transformed metric: 



M = -xM^<i> + Six); 
ax 



(24) 



52^ 

dxl 



(28) 



where ^ is given by H21I) and we are once again con- 
sidering periodic perturbations. As might be expected, 
the potential C is more complicated than for the neutral 
black holes: 

C = ^ [1 - NA^^B^ - 2Ar^ - iSar^ - Ylar^A 
3 

3or^ a \dr 



d(t) 
dr 



A-'B-^ 



60 a \ dr 



+Sar^^^^A-'B-' 
dr 



lAr^^A-'B-' 
3 dr 



(29) 



where S is an arbitrary function of x. With these ex- 
pressions, the (tt) component of the linearized Einstein 
equation gives an equation of the form 

XF + y^ + gsix)+n^^^ = 0; (25) 

ox 

where X, y, Q and Ti are functions of x only. Using 
the equilibrium field equations, it can be shown that X 
vanishes, and 



y = -2— + x( — 

ox \ Ox 



(26) 



Furthermore, it can be shown that 3^ = by writing the 
right-hand-side of H26|l in terms of the original variables 
and radial coordinate r. Equation (|25|l can then be inte- 
grated to give 



S = K exp — 



-dx 



K 



with A and B as before H19(l . As in the case of neutral 
black holes, the Schrodinger equation (|28|l with poten- 
tial could have been derived directly without using 
the conformal transformation l|12|l. However, our method 
has made the computations required significantly more 
tractable. 

The potential C H29|l reduces, in the limit a — > —A/36, 
to the potential U for the neutral black holes. This is 
as expected, since in this limit the charge Q of the black 
hole tends to zero. 

It is again straightforward to show that d(j)/dr and 
B are all regular and non-zero everywhere between the 
event and cosmological horizons. However, unlike in the 
neutral case, for the charged MTZ black holes it is no 
longer necessarily the case that A vanishes for some r 
lying between r_|_ and Putting (p = \/6 at r = ro 

in (IHI and using the inequalities M < |; r+ < tq < ?'++ 
and < 3^ < 1, together with 0, gives the condi- 
tion that A vanishes somewhere between the event and 
cosmological horizons if and only if 



where K is an arbitrary constant. Now, using the bound- 
ary conditions that F = $ = at the event horizon, it 
must be the case that K — and so S{x) — identically. 
Using this result, the (xx) component of the linearized 
Einstein equations reduces to a relationship between the 
electromagnetic field and scalar field perturbations: 



F 



dF dF 
dx dx 



F = xF- 



9$ 9$ 

dx dx 



(27) 



M < (-12a)- 




(30) 



When Q = we have from (0 that a = and 
becomes M < j, which is always satisfied, retrieving 
the result of section IHII Similarly, in the limit A 
0, the inequality (|30|) is trivially satisfied for all M, in 
accordance with the results for the BBMB black hole @| . 
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FIG. 2: The region 3?, described by Eq. Inside K, 

the potential C has a pole (see Fig. , and the corresponding 
black holes are unstable. Outside K, the potential C is regular 
(see Fig. I^J, and numerical methods are necessary to show 
that the black holes are unstable. 

If equation (|30|l is satisfied, then the second term in 
the potential C (^5)) has a double pole at r = ro (see Fig. 
nfor a typical potential in this case). Although there are 
other terms in C which have a single pole at r = ro, these 
will be sub-leading compared to the double pole. In this 
case the behaviour of C near r = rp is the same as for the 
potential U (tTH|) : 

C = -\ [N{ro)f (r - ro)"' + O (r - ro)"' . 

As in the analysis of the neutral black holes, in this case 
standard results in quantum mechanics Q allow us to de- 
duce that there are an infinite number of unstable modes. 

It remains therefore to understand whether the in- 
equality H3U|I is satisfied by any or all of the charged MTZ 
black hole solutions. We express this inequality in terms 
of the fractions of the maximum possible values of M 
(1/4:) and Q {l^jy/^ for any particular value of A, and 
denote these model parameters by m and q respectively. 
The region Sft of the (m, q) phase space for which we have 
shown the solution to be unstable is then 

5R=|m(A),g(A),A:m<-iil^^^ll-|. (31) 

[ (l + (l-g2),)2j 

The region is shown in figure [3 where it can be seen 
that it makes up the vast majority of the available param- 
eter space. Only black holes carrying a very high charge, 
or those nearly as large as their universe (so that the 
event and cosmological horizons nearly coincide) , are pa- 
rameterized by variables (m, q) lying outside the region 
3?. While it could be argued that such black holes are 
unphysical, nevertheless for completeness we investigate 
the phase space outside 3?. 

A typical potential for a point outside 3? is shown in 
Fig. 13 The potential is regular everywhere between 



Potential for a black hole configuration outside 9t 




Radius 



FIG. 3: Form of the potential C for charged MTZ black hole 
configurations given by parameters (m, q) lying outside the 
region K. The potential is regular everywhere between the 
event and cosmological horizons. It is positive close to the 
event horizon, and then becomes negative. As the values of 
(m, q) approach the boundary of the region K, the positive 
peak close to the event horizon becomes larger and larger. 



the event and cosmological horizons, where it vanishes. 
There is a region close to the event horizon where the 
potential is positive; elsewhere it is negative. For black 
holes described by parameters (m, q) lying close to the 
boundary of the region 3?, the positive peak close to the 
event horizon becomes very large. 

We show in the appendix that potentials of this form 
with 

/oo 
Cda;, <0; (32) 
-oo 

must have at least one bound state, and therefore the 
corresponding black holes will be unstable. We have cal- 
culated numerically the integral on the left-hand-side of 
(I32|l ■ and some typical results are shown in Fig. 01 We 
find that for most configurations outside 3?, the integral 
(|32|l is indeed negative, so that those black holes are un- 
stable. 

However, for black hole configurations very close to the 
boundary of 3?, we find that the integral becomes posi- 
tive (due to the large positive peak in the potential close 
to the event horizon). For such black holes their stabil- 
ity has to be investigated numerically, each configuration 
being considered separately. A simple way to check the 
existence of negative eigenvalues of (|28|) . without actually 
finding the values of these eigenvalues, is to apply conti- 
nuity arguments of the type outlined in Ref. [iq . If, for 
some (T^ < (which is not necessarily an eigenvalue of 
(|28|l '). the solution of (|28|l satisfying the boundary con- 
ditions at, say, a;* — > — oo (r = r_|_) has one zero before 
possibly diverging as x» ^ oo (r ^ then there 

is at least one negative eigenvalue of (|28|l . We studied 
various black holes corresponding to points (m, q) lying 
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FIG. 4: The integral of the potential C at points lying outside 
the region 5R. In this case the charge is fixed to be 0.95 of the 
maximum charge, and the values of the mass vary from the 
value on the boundary of 5R up to the maximum. 



outside the region and such that the integral on the 
left-hand-side of (|32|) is positive, and found, in each case, 
a value of fi'^ < satisfying this condition. Therefore, for 
all the cases we studied, there are negative eigenvalues of 
(|28|l . showing instability. 



V. CONCLUSIONS 

We have investigated the stability of the de Sitter black 
hole solutions of Einstein's equations due to Martinez et 
al 0. The model contains a conformally coupled scalar 
field with a quartic self-interaction potential and also an 
electromagnetic field in the case that the black holes are 
charged. We have shown by analytic methods that all the 
neutral black holes are unstable, as are the charged black 
holes in the vast majority of the phase space. Outside this 
region of phase space, we have used numerical methods 
to show that the black holes are unstable. 

The analysis differs from that for the analogous, 
asymptotically flat BBMB black holes, where the 
analytic approach suffices for both the neutral and all 
the charged configurations. However, the conclusions 
are the same, that the neutral and charged MTZ black 
holes, like the BBMB black holes are unstable. We are 
therefore able to complete the following table, describing 
the behaviour of conformally coupled scalar field black 
hole hair (cf. that in Ref. Q): 



A = Unstable hair 
A > Unstable hair 
A < Stable hair 
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APPENDIX: INSTABILITY IF THE INTEGRAL 
OF THE POTENTIAL IS NEGATIVE 

In this appendix we shall prove a result needed in the 
stability analysis of the charged MTZ solutions in section 
IIVI namely that if 



C dx^ < 0; 



(A.l) 



then there is at least one negative eigenvalue of H28|l , and 
the corresponding black hole configuration is unstable. 

The starting point is a result from Ref. if 
there exists a twice differentiable function / such that 
/(-oo) = /(oo) = 0, and 



-7^ +cr 



dxt < 0; 



(A.2) 



then there is at least one negative eigenvalue of H28|l . 

In this appendix we show that the condition (|A.2|I is 
satisfied if (|A.1(I holds. The proof is straightforward, so 
we shall just briefiy outline the key steps. 

We define a sequence of functions fk{x*) : fc € N, such 
that each fk{x>,) tends to zero as |a;*| oo, following 
M: 



where 



and 



Ziu) ~ Z{—u); 

Z{u) ^1 for ue [0,p]; 

Z{u) = for u> p+1; 



-q<Z'{u)<0 for ue[p,p+l]; 



(A.3) 



where p and q are positive constants. Therefore fk{x^,) is 
an even function stretched horizontally for increasing k. 
Now is non-zero only on the intervals 

h ^ [-k{p+l),~kp]; h = [kp,k{p+l)]- 

and decreases linearly with k for corresponding points as 
fk is scaled. Hence correspondingly decreases as 
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■p- , whereas the widths of the intervals Ii and I2 increase 
hnearly with k. It follows that 




decreases proportionally as as increases, hence the 
first term in ljA.2|) can be made arbitrarily small for suf- 
ficiently large k. 

We know that for black hole configurations correspond- 
ing to points lying outside the region (defined in section 
II Vp . the potential C is continuous and tends to zero as 
X* —^ ±00, or as r — > r+ or r++. From this and the 
definition of fk{x*) HA.3|) we deduce that 

/oo poo 
Cdx^ — i {fk{x*)) Cdx^, 
-oo J ~oo 

/—pk 1*00 
\C \ dx^ + \C \ dx-t. 
-00 J pk 



In addition, the potential C approaches zero linearly 
with respect to r very near and r_|__|_, so equation 
(fTTjl gives us that for large \x^\, the potential \C\ de- 
creases to zero exponentially with Hence, both 

S-^ |C| dx^ and \C\ dx^ decrease exponentially with 
increasing k. Therefore, /^Cc?x* will be arbitrarily 
close to C/^ dx* for sufficiently large k. This gives 
us our result, that if Cdx* is negative we can pick a 
sufficiently large k such that equation ljA.2|) holds with 
f = fk, and we have shown the solution to be unstable. 
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